Spin-orbit precession damping in transition metal ferromagnets 
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We provide a simple explanation, based on an effective field, for the precession damping rate due 
to the spin-orbit interaction. Previous effective field treatments of spin-orbit damping include only 
variations of the state energies with respect to the magnetization direction, an effect referred to 
as the breathing Fermi surface. Treating the interaction of the rotating spins with the orbits as 
a perturbation, we include also changes in the state populations in the effective field. In order to 
investigate the quantitative differences between the damping rates of iron, cobalt, and nickel, we 
compute the dependence of the damping rate on the density of states and the spin-orbit parameter. 
There is a strong correlation between the density of states and the damping rate. The intraband 
terms of the damping rate depend on the spin-orbit parameter cubed while the interband terms are 
proportional to the spin-orbit parameter squared. However, the spectrum of band gaps is also an 
important quantity and does not appear to depend in a simple way on material parameters. 



I. INTRODUCTION 

Magnetic memory devices are useful if they can be re- 
liably switched between two stable states. The fidelity of 
this switching process depends sensitively on the damp- 
ing rate of the system. Despite decades of research and 
the relentless industrial push toward smaller and faster 
devices, many questions about the damping process re- 
main unanswered, particularly for metallic ferromagnets. 
Recent experimental efforts have investigated the extent 
to which the damping rate of NiFe alloys can be tuned 
through doping, particularly with the addition of rare 
earth [l[ and transition metal elements While these 
investigations found a general trend suggesting damp- 
ing increases with increasing spin-orbit coupling of the 
dopant, the details behind this effect remained elusive. 
To aid this effort, this article provides a simple descrip- 
tion of the damping process and investigates how some 
material properties affect the damping rate. 

Precession damping in metallic ferromagnets results 
predominantly from a combined effort of spin-orbit cou- 
pling and electron- lattice scattering 0, 0]. The role of 
lattice scattering was studied in early experimental work 
through the temperature dependence of damping rates 
[H, @ . Measurement of damping rates versus tempera- 
ture revealed two primary contributions to damping, an 
expected part that increased with temperature, and an 
unexpected part that decreased with temperature. In 
cobalt these two opposing contributions combine to pro- 
duce a minimum damping rate near 100 K, for nickel the 
increasing term is weaker leading to a temperature inde- 
pendent damping rate above 300 K, while for iron the 
damping rate becomes independent of temperature be- 
low room temperature. Heinrich et al. later noted that 
the temperature dependence of the increasing and de- 
creasing contributions matched that of the resistivity and 
conductivity, respectively @, [1] , and so dubbed the two 
contributions conductivity-like for the decreasing piece 
and resistivity- like for the increasing part. 

Among the many theories on intrinsic precession 



damping 0, i ES El El El El MM, El El P , Kam- 
bersky's torque-correlation model [3J is unique in qualita- 
tively matching the observed non-monotonic temperature 
dependence that we just described. We recently evalu- 
ated this model for iron, cobalt, and nickel, and showed 
that it accurately predicts the precession damping rates 
of these systems [4j . While this model succeeds in captur- 
ing the important physical effects involved in precession 
damping, it does not easily identify the important physi- 
cal processes or give insight into how one might alter the 
damping rate through sample manipulation. In section II 
we briefly describe the torque-correlation model. In or- 
der to provide a more tangible explanation of precession 
damping we rederive the damping rate from an effective 
field approach in section III. This discussion is followed 
in section IV by a quantitative analysis of the effect on 
the damping rate of tuning the density of states and the 
spin-orbit parameter. 



II. TORQUE-CORRELATION MODEL 

Kambersky's theory describes damping in terms of the 
spin-orbit torque correlation function, finding a damping 
rate of 



x / de 1 A nk (e 1 )A mk (e 1 )i](e 1 ) . (1) 



The gyromagnetic ratio is 7 = gnofiB/h, 9 is the Lande g 
factor, no is the permeability of space, n and m are band 
indices, and k is the electron wavevector. The matrix ele- 
ments |r~ n (/c)| 2 describe a torque between the spin and 
orbital moments that arises as the spins precess. 77(e) 
is the derivative of the Fermi function —df/de, which 
is a positive distribution peaked about the Fermi level 
that restricts scattering events to the neighborhood of 
the Fermi surface. The electron spectral functions A n k(e) 
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are Lorentzians in energy space centered at band energies 
with widths determined by the scattering rate. They phe- 
nomenologically account for electron-lattice scattering. 

Equation ([!} includes two processes: the decay of 
magnons into electron-hole pairs and the scattering of 
the electrons and holes with the lattice. This expres- 
sion is similar in structure to sp-d models that have 
proven successful in describing dissipation in semicon- 
ductors • However, the physics of the magnon decay 
process is very different. In the present case, there is no 
distinction between sp and d electrons. The spin-orbit 
torque annihilates a uniform mode magnon and gener- 
ates an electron-hole pair. The electron-hole pair is then 
collapsed through lattice scattering. The electron and 
hole are dressed through lattice interactions and are best 
thought of as a single quasiparticle with indeterminant 
energy and a lifetime given by the electron-lattice scat- 
tering time. The dressed electron and hole can occupy 
the same band (to = n), which we call an intraband 
transition, or two different bands (to ^ n), an inter- 
band transition. For intraband transitions, the integra- 
tion over the spectral functions is proportional to the 
scattering time, just like the conductivity. For intcrband 
transitions, the intregration over the spectral functions 
is roughly inversely proportional to the scattering time, 
as is the resistivity. Therefore, the intraband terms in 
Eq. (fTJ give the conductivity-like contributions to damp- 
ing that decrease with temperature while the interband 
terms yield the resistivity-like contributions that increase 
with temperature. 



III. EFFECTIVE FIELD DERIVATION 

An effective field for the magnetization dynamics is de- 
fined as the variation of the electronic energy with respect 
to the magnetization direction /xoH cff = —dE/dNL. The 
magnitude of the magnetization M is considered con- 
stant within the Landau-Lifshitz formulation, only the 
direction M of the magnetization changes. The total 
electronic energy of the system can be approximated by 
E = ^2 nk Pnktnk, which is a summation over the single 
electron energies e n k weighted by the state occupancies 
p n k ■ If the state occupancies are held at their equilbrium 
values, the resulting effective field is equivalent to that of 
the magnetocrystalline anisotropy [21(, which describes 
reversible processes. If however, the state occupancies 
are allowed to deviate from the equilibrium populations 
in response to the oscillating perturbation, an irreversible 
contribution also arrises, which we show produces the 
damping in Eq. JTj). 

As the magnetization precesses the energies of the 
states change through variations in the spin-orbit con- 
tribution and transitions between states occur. These 
two effects, the changing energies of the states and the 
transitions between states, produce a contribution to the 
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The first term in the brackets describes the variation in 
the spin-orbit energies of the states as the magnetization 
direction changes. This effect, which has been discussed 
and evaluated before [H, El, is generally referred 
to as the breathing Fermi surface model. The spin-orbit 
torque does not cause transitions between states in this 
picture, but does cause the Fermi surface to swell and 
contract as the magnetization precesses. We will show 
that this portion of the effective field gives the intraband 
terms of Eq. |T]). The second term in the brackets has pre- 
viously been neglected in effective field treatments, but 
accounts for changes in the system energy due to tran- 
sitions between states. This term does not change the 
energies of the states, but does create electron-hole pairs 
by exciting electrons from lower bands to higher bands. 
This process can be pictured as a bubbling of individ- 
ual electrons on the Fermi surface. We will demonstrate 
that this portion of the effective field gives the interband 
terms of Eq. (JTJ). 



A. Intraband terms 

The first term in the effective field Eq. @ accounts for 
the effects of the breathing Fermi surface (bfs) model. 
Since this model has previously been discussed in detail 
[Tol . IT3I l22l | we will give only a very brief review of it here, 
focusing instead on connecting it to the intraband terms 
of Eq. 0. 

As the magnetization precesses the spin-orbit energy 
of each state changes. Some occupied states originally 
just below the Fermi level get pushed above the Fermi 
level and simultaneously some unoccupied state origi- 
nally above the Fermi level may be pushed below it. 
This process takes the system, which was originally in 
the ground state, and drives it out of equilibrium into an 
excited state creating electron-hole pairs in the absence 
of any scattering events. Scattering, which occurs with a 
rate given by the inverse of the relaxation time r, brings 
the system to a new equilibrium. The relaxation time 
approximation determines how far from equilibrium the 
system can get. 



Pnk — fnk 



_ dfnk 

dt 



(3) 



The occupancy p n k of each state ipnk deviates from its 
equilibrium value f n k by an amount proportional to the 
scattering time. How quickly the system damps depends 
on the magnitude of this deviation. 

The rate of change of the equilibrium distribution 
dfnk/dt depends on how much the distribution changes 
as the energy of the state changes df n k/de n f., how much 
the state energy changes as the precession angle changes 
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(a) 



FIG. 1: Schematic description of precession geometry. Within 
the breathing Fermi surface model (a) the damping rate is cal- 
culated as the magnetization passes through a specific point in 
a given direction. The torque correlation model (b) gives the 
damping rate for precessing about a given direction. Dashed 
curves indicate the precession trajectory. 



de n k/dM, and how quickly the spin direction is precess- 
ing dM/dt. These can be combined with a chain rule 



dfnk df nk de nk dM 



dt 



dtnk dM dt 



(4) 



Combining this result with the relaxation time approxi- 
mation Eq. ([3]) and substituting these state occupancies 
into the first term of the effective field in Eq. ([2]) gives 
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H™g is a contribution to the magnetocrystalline 
anisotropy field and H^j s mp the damping field from the 
breathing Fermi surface model. When we compare this 
damping field to the damping field postulated by the 
Landau-Lifshitz-Gilbert equation 



XT-damp 



A dM 
'j 2 M dt 



(8) 



we find that the damping rate is 



model Eq. (P) gives the damping rate when the magne- 
tization is undergoing small angle precession about the 
z direction (see Fig.(jT|b)). When z is a high symmetry 
direction the change in the magnetization will stay in the 
x-y plane. In each scenario - rotating M through z in the 
breathing Fermi surface model and rotating M about z 
in the torque correlation model - dM is confined to the 
x-y plane. Therefore, rotating through z and rotating 
about z are equivalent in the small angle limit when z 
is a high symmetry direction. With this observation we 
now show that the intraband contributions of the torque 
correlation model are equivalent to the breathing Fermi 
surface result under these conditions. 

The only energy that changes as the magnetization 
rotates is the spin-orbit energy 7i so . As the spin of the 
state \nk) rotates about the d direction by angle d its 
spin-orbit energy is given by 



e(i?) = (nk\ 



ik) 



(10) 



where § = Taking the derivative of this energy with 
respect to $ in the limit that goes to zero shows that 
the energy derivatives are 



de 

— = i(nk\{a ■ ft, H so ]\nk) 



(11) 



Figure ([I]) shows that the derivative de/dd is identical to 
de/dM and that when M = z the rotation direction lies 
in the x — y plane. The two components of the transverse 
torque operator T x and T v can be obtained (up to factors 
of i) by setting d equal to x or y, respectively. From this 
observation we find 



|(nfc|r-|nfc)| 2 = 



0c_ 

Ox 



(12) 



When the magnetization direction z is pointed along a 
high symmetry direction the transverse directions x and 
y are equivalent and | r | 2 = 2(de/dM) 2 . 

Substituting the torque matrix elements for the energy 
derivatives in Eq.© gives a damping rate of 



Abfs = t — > r){e n k) 



Mo 



dM 



(9) 



\bfs 



T7 

2^o 



(13) 



n k 



As in Eq. ([T]), 77(e) is the negative derivative of the Fermi 
function and is a positive distribution peaked about the 
Fermi energy. 

As described in Fig. ([IK), the result of the breathing 
Fermi surface model Eq. © describes the damping rate 
of a material as the magnetization rotates through a par- 
ticular point z about a given axis t9. When M is instan- 
taneously aligned with £ the direction of the change in 
the magnetization dM will be perpendicular to z, in the 
x-y plane. On the other hand, the torque correlation 



For the intraband terms in Eq. |T]) the integration over 
the spectral functions reduces to rr](e n k)/2Trh so we find 



Hits 



Mo ^ 



d 3 k 



K(k)\ 



(2tt) 3 

deiA„ fc (ei)A„ fe (ei)?7(ei) , 
which matches the intraband terms of Eq. |T|). 



(14) 
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B. Interband terms 

As the magnetization precesses, the spins rotate and 
the spin-orbit energy changes. This variation acts as a 
time dependent perturbation 

V(t) = e^^Hsoe'^^ - H so (0) rj i[a ■ <p(t),H so ] . 

This approximation results from linearizing the expo- 
nents, which is appropriate in the small angle limit. 
The time dependence of the spin direction is (p(t) — 
cos cot x + shxcoty, up to a phase factor. This perturba- 
tion causes band transitions between the states ip nk and 
ijjmk- The initial and final states have the same wavevec- 
tor because these transitions are caused by the uniform 
precession, which has a wavevector of zero. The transi- 
tion rate between states due to this perturbation is 

2-7T | ,2 

W mn (k) = — T mn (fc) S(e mk - e nk - Hoj) . (16) 
h 1 1 

The variations of the occupancies of the states with 
respect to the magnetization direction are given by the 
master equation 



dt 



^ W mn (k)[p m k - Pnk] ■ (17) 



The second term in the effective field Eq. (|2J) contains 
the factor dp nk /dM which is (dp nk /dt)/(dip/dt) where 
dtp/dt — io. Inserting these expressions into the second 
term in the effective field and rearranging the sums gives 



H 



eff_ 1 \ - \ - W mn (k) t dM 
2flQ M^^ to 2 [Pnk Pmk ^ tmk €nk > dt 

nk m^n 

(18) 

Comparing this result to the effective field predicted by 
the Landau-Lifshitz-Gilbert equation (jHJ) we find a damp- 
ing rate of 

\ 7 V~~* \ ^ jjrr /i \ [Pnk ~ Pmk] [ e mk ~ e nk\ f-ir\\ 

A = — > > W mn {k) . 19) 

nk ra^n 

The finite lifetime of the states is introduced with the 
spectral functions 



A 



2/i 

x W mn (k) 
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[/(ci)-/(e 2 )] [e 2 -ei] 



fi/jo 



(20) 



Inserting the transition rate Eq. (I16|l . integrating over €2, 
and taking the limit that to goes to zero leaves 



A 



4EE 



/to 



d 3 k 



n m^n ' 

x / de 1 A nk (e 1 )A mk (e 1 )r](e 1 ) , 



(21) 



which are the interband terms of Eq. ([T]) . 

In this derivation of the bubbling Fermi surface contri- 
bution to the damping we have ignored an additional, re- 
versible term that contributes to the magnetocrystalline 
anisotropy. This contribution arises from changes in the 
equilibrium state occupancies as the magnetization di- 
rection changes. This contribution to the magnetocrys- 
talline anisotropy is localized to the Fermi surface while 
the contribution dicussed in IIIA is spread over all of the 
occupied levels. 



IV. TUNING THE DAMPING RATE 

We have previously demonstrated that the mcchansim 
of thetorque correlation model Eq. (fTJ) accounts for the 
majority of the precession damping rates of the transi- 
tion metals iron, cobalt, and nickel In the present 
work we have so far shown that this expression for the 
damping rate can be described simply within an effective 
field picture. We now investigate the degree to which 
the damping rate may be modified by adjusting certain 
material parameters. Inspection of Eq. fT]) reveals that 
the damping rate depends on the convolution of two fac- 
tors: the torque matrix elements and the integral over the 
spectral functions. We separate the quantitative analy- 
sis of the damping rates into their dependencies on these 
two factors, beginning with the spectral weight. 

The calculations for the damping rate of Eq. (TT|) dis- 
cussed below are performed using the linear augmented 
planewave method in the local spin density approxima- 
tion. The details of the computational technique may be 
found in H, |2ll. and the included references. 



A. Spectral overlap 

For the intraband terms, the integral over the spec- 
tral functions is essentially proportional to the density of 
states at the Fermi level. Therefore, it appears reason- 
able to suspect that the intraband contribution to the 
damping rate of a given material should be roughly pro- 
portional to the density of states of that material at the 
Fermi level. To test this claim numerically, we artificially 
varied the Fermi level of the metals within the d-bands 
and calculated the intraband damping rate as a function 
of the Fermi level. The results of these calculations are 
superimposed on the calculated densities of states of the 
materials in Fig. O The correlation between the damp- 
ing rates and the densities of states, while not exact, is 
certainly strong, indicating that increasing the density 
of states of a system at the Fermi level will generally 
increase the intraband contribution to damping. 

The dependence of the interband terms on the spectral 
overlap is more complicated than that of the intraband 
terms. The spectral overlap depends on the energy dif- 
ferences e m — e„, which can vary significantly between 
bands and over k-points. When the scattering rate H/r 
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FIG. 2: Intraband damping rate versus Fermi level superim- 
posed upon density of states. A strong correlation between 
the intraband damping rate versus Fermi level (•) and the 
density of states (solid curves) is observed. Vertical black 
lines indicates true Fermi energy calculated by density func- 
tional theory. 



is much less than these energy gaps the interband terms 
are proportional to the scattering rate. However, this 
proportionality only holds at low scattering rates when 
the interband contribution is much less than the intra- 
band contribution. The proportionality breaks down at 
higher scattering rates when hjr becomes comparable to 
the band gaps. After this point the damping rate gradu- 
ally plateaus with respect to the scattering rate. Unfor- 
tunately, this complicated functional dependence of the 
spectral overlap on the scattering rate makes it difficult 
to obtain a simple description of the effect of the spec- 
tral overlap on the interband damping rate in terms of 
material parameters. 



B. Torque matrix elements 

The damping rate also depends on the square of the 
torque matrix elements. A goal of doping is often to mod- 
ify the effective spin-orbit coupling of a sample. While 
doping does more than this, such as intoducing strong lo- 
cal scattering centers, it is nevertheless useful to estimate 
the dependence of the matrix elements on the spin-orbit 
parameter £. We begin with pure spin states ip^ and 
treat the spin-orbit interaction V = £V as a perturba- 
tion. The states can be expanded in powers of £ as 

^=< + ^ + £ 2 ^ + ---. (22) 

The superscripts refer to the unperturbed wavefunction 
(0) and the additions (i) due to the perturbation to the 
ith order while the subscript n is the band index, which 
includes the spin direction, up or down. Since the torque 
operator also contains a factor of the spin-orbit parame- 
ter the matrix elements have terms in every order of £ be- 
ginning with the first order. Therefore, the squared ma- 
trix elements have contributions of order £ 2 and higher. 

To determine the importance of these terms we arti- 
ficially tune the spin-orbit interaction from zero to full 
strength, calculating the damping rate over this range. 
We then fit the intraband and interband damping rates 
separately to polynomials. In each material, this fitting 
showed that for the intraband terms the £ dependence of 
the damping rate was primarily third order, with smaller 
contributions from the second and fourth order terms. 
Restricting the fit to only the third order term produced 
a very reasonable result, shown in Fig. ([3]). For the in- 
terband terms, polynomial fitting was dominated by the 
second order term, with all other powers contributing 
only negligibly. The second order fit is shown in Fig. ([3]). 

To understand the difference in the £ dependence of 
the intraband and interband contributions it is useful to 
define the torque operator 

T- =£,{ro z -ta-). (23) 

The torque operator lowers the angular momentum of 
the state it acts on. This can be accomplished either by 
lowering the spin momentum l z o~ , a spin flip, or low- 
ering the orbital momentum l~ o z , an orbital excitation. 
Therefore, both the intraband and interband contribu- 
tions each have two sub-mechanism: spin flips and orbital 
excitations. 

The second order terms for the intraband case are 
£ 2 1 (V'n I a ~~ ^ ° z ) I ) 1 2 ■ Since the unperturbed states 
are pure spin states the spin flip part £ z a~ of the 
torque returns zero. Therefore, only the orbital excita- 
tions exist to lowest order in £, reducing the strength 
of the second order term in the intraband case. How- 
ever, the interband terms contain matrix elements be- 
tween several states, some with the same spin direction, 
but others with opposite spin direction. Therefore, both 
spin flips and orbital excitations contribute in second or- 
der to the interband contribution. 
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FIG. 3: £ dependence of intraband and interband damping 
rates. Damping rates were calculated for a range of spin-orbit 
interaction strengths between off (£ = 0) and full strength 
(£ = 1). £ 2 fits were made to the interband damping rates 
(left axes and ■< symbols) and £ 3 fits to the intraband rates 
(right axes and ► symbols). 



V. CONCLUSIONS 

The breathing Fermi surface model has provided a sim- 
ple and understandable effective field explanation of pre- 
cession damping in metallic ferromagnets. However, it is 
only applicable to very pure systems at low temperatures. 
On the other hand, the torque correlation model accu- 
rately predicts damping rates of systems with imperfec- 
tions from low temperatures to above room temperature. 
The shortcoming of the torque correlation model is that 
it does not illuminate the phyiscal mechanisms responsi- 
ble for damping. We have pointed out that the breath- 
ing Fermi surface model accounts for only one of the two 
terms in the effective field. By constructing an effective 
field with the previously studied breathing Fermi surface 
contribution and also the new bubbling effect we have 



shown that this simpler picture may be mapped onto the 
torque correlation model such that the breathing terms 
match the intraband contribution and the bubbling terms 
match the interband contribution. 

Since there is considerable interest in understanding 
how to manipulate the damping rates of materials we 
investigated the dependence of the intraband and inter- 
band damping rates on both the spectral overlap inte- 
gral and the torque matrix elements. For the intraband 
terms, the spectral overlap is proportional to the density 
of states and we found a strong correlation between the 
intraband damping rate and the density of states of the 
material. The interband case is significantly complicated 
by the range of band gaps present in materials. No sim- 
ple relation was found between the strength or scattering 
rate dependence of the interband terms and common ma- 
terial parameters. The importance of the torque matrix 
elements to the damping rates was characterized through 
their dependence on the spin-orbit parameter. The intra- 
band damping rates were found to vary as the spin-orbit 
parameter cubed while the interband damping rates went 
as the spin-orbit parameter squared. This difference was 
explained by noting that the torque operator changes the 
angular momentum of states either through spin flips, or 
by changing their orbital angular momentum. Spin-flip 
excitations do not occur to second order in £ for the in- 
traband terms, but do contribute at second order for the 
interband terms. 

It is desirable to understand the relative differences in 
damping rates amoung various materials, such as why 
the damping rate for nickel is higher than that for cobalt 
and iron. We have shown that the relative damping rates 
of these materials depend in part on the differences of 
their densities of states and spin-orbit coupling strengths. 
However, they also depend in an intricate way on the en- 
ergy gap spectra of each metal. For the interband terms 
the dependence on the gap spectrum enters through the 
spectral overlap integral. For the intraband terms the 
energy gaps appear in the denominators of the matrix 
elements. Therefore, states with very small splittings 
can dominate the k-space convolution. The abundance 
of such states in nickel appears to contribute to the larger 
damping rate in this material [23| . 

Doping is a common technique for modifying damp- 
ing rates. Doping has a number of consequences on a 
sample and these effects vary with the method of dop- 
ing. Dopants can increase the electron-lattice scattering 
rate, introduce magnetic inhomogeneities that act as lo- 
cal scattering centers, alter the density of states, and 
change the effective spin-orbit parameter. We have in- 
vestigated the consequences of modifying the densities of 
states and spin-orbit parameter on the damping rate, and 
previously demonstrated the scattering rate dependence 
of the damping rate; however, it is not clear what new 
damping mechanisms arise when rare-earth elements are 
added to a transition metal host. 

This work was supported in part by the Office of Naval 
Research through grant N00014-03-1-0692 and through 
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